(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 7.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 43758, 1095] NotebookOptionsPosition[ 41222, 1001] NotebookOutlinePosition[ 41678, 1020] CellTagsIndexPosition[ 41635, 1017] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "\:041f\:0435\:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\:0430\:044f", " ", "\:0444\:0443\:043d\:043a\:0446\:0438\:044f", " ", "\:0441\:043a\:043e\:0440\:0435\:043a\:0442\:0438\:0440\:043e\:0432\:0430\ \:043d\:043d\:043e\:0439", " ", "\:0441\:0438\:0441\:0442\:0435\:043c\:044b", " ", "\:0434\:0435\:043c\:043f\:0444\:0435\:0440\:0430", " ", "\:0440\:044b\:0441\:043a\:0430\:043d\:0438\:044f", " ", RowBox[{"Wy", "/", "\[Delta]\:043d"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"w", "=", FractionBox[ RowBox[{"1", "*", "7.56", " ", RowBox[{"(", RowBox[{"0.47", "\[InvisibleSpace]", "+", "s"}], ")"}], " ", RowBox[{"(", RowBox[{"0.47", "\[InvisibleSpace]", "+", "s"}], ")"}]}], RowBox[{ RowBox[{"(", RowBox[{"0.47", "\[InvisibleSpace]", "+", "s"}], ")"}], " ", RowBox[{"(", RowBox[{"41.01", "\[InvisibleSpace]", "+", RowBox[{"12.80", " ", "s"}], "+", SuperscriptBox["s", "2"]}], ")"}]}]]}], "\[IndentingNewLine]", RowBox[{"Wt", "=", FractionBox[ RowBox[{"1", "*", "7.56", " ", RowBox[{"(", RowBox[{"0.47", "\[InvisibleSpace]", "+", "s"}], ")"}], " ", RowBox[{"(", RowBox[{"0.47", "\[InvisibleSpace]", "+", "s"}], ")"}]}], RowBox[{ RowBox[{"(", RowBox[{"0.47", "\[InvisibleSpace]", "+", "s"}], ")"}], " ", RowBox[{"(", RowBox[{"41.01", "\[InvisibleSpace]", "+", RowBox[{"12.80", " ", "s"}], "+", SuperscriptBox["s", "2"]}], ")"}]}]]}]}]}]], "Input"], Cell[BoxData[ FractionBox[ RowBox[{"7.56`", " ", RowBox[{"(", RowBox[{"0.47`", "\[InvisibleSpace]", "+", "s"}], ")"}]}], RowBox[{"41.01`", "\[InvisibleSpace]", "+", RowBox[{"12.8`", " ", "s"}], "+", SuperscriptBox["s", "2"]}]]], "Output"], Cell[BoxData[ FractionBox[ RowBox[{"7.56`", " ", RowBox[{"(", RowBox[{"0.47`", "\[InvisibleSpace]", "+", "s"}], ")"}]}], RowBox[{"41.01`", "\[InvisibleSpace]", "+", RowBox[{"12.8`", " ", "s"}], "+", SuperscriptBox["s", "2"]}]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "\:041f\:0435\:0440\:0435\:0439\:0434\:0435\:043c", " ", "\:043e\:0442", " ", "s", " ", "\:043a", " ", "\:043a\:043e\:043c\:043f\:043b\:0435\:043a\:0441\:043d\:043e\:0439", " ", "\:0444\:043e\:0440\:043c\:0435"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"wk", "=", RowBox[{"w", "/.", RowBox[{"s", "\[Rule]", RowBox[{"I", "*", "omega"}]}]}]}]}]], "Input", CellChangeTimes->{{3.5160879222637167`*^9, 3.5160879667549953`*^9}, { 3.516088011917074*^9, 3.5160880449579325`*^9}}], Cell[BoxData[ FractionBox[ RowBox[{"7.56`", " ", RowBox[{"(", RowBox[{"0.47`", "\[InvisibleSpace]", "+", RowBox[{"\[ImaginaryI]", " ", "omega"}]}], ")"}]}], RowBox[{"41.01`", "\[InvisibleSpace]", "+", RowBox[{"12.8`", " ", "\[ImaginaryI]", " ", "omega"}], "-", SuperscriptBox["omega", "2"]}]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"f", "=", RowBox[{ RowBox[{"(", RowBox[{"Abs", "[", "wk", "]"}], ")"}], "^", "2"}]}]], "Input", CellChangeTimes->{{3.5160881542985244`*^9, 3.5160881791025677`*^9}, { 3.5160882708931293`*^9, 3.5160882833419514`*^9}}], Cell[BoxData[ RowBox[{"57.1536`", " ", SuperscriptBox[ RowBox[{"Abs", "[", FractionBox[ RowBox[{"0.47`", "\[InvisibleSpace]", "+", RowBox[{"\[ImaginaryI]", " ", "omega"}]}], RowBox[{"41.01`", "\[InvisibleSpace]", "+", RowBox[{"12.8`", " ", "\[ImaginaryI]", " ", "omega"}], "-", SuperscriptBox["omega", "2"]}]], "]"}], "2"]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"fvhod", "=", "1", " ", RowBox[{"(*", RowBox[{ "\:0432\:0445\:043e\:0434\:043d\:043e\:0439", " ", "\:0441\:043f\:0435\:043a\:0442\:0440"}], "*)"}]}]], "Input", CellChangeTimes->{{3.5160882915163655`*^9, 3.5160882979903765`*^9}, { 3.516089737672906*^9, 3.5160897458785205`*^9}, {3.516343688283702*^9, 3.516343691169707*^9}}], Cell[BoxData["1"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"fvihod", "=", RowBox[{"f", "*", "fvhod"}]}]], "Input", CellChangeTimes->{{3.5160883009387817`*^9, 3.5160883169444103`*^9}, { 3.516343695818515*^9, 3.5163436974097176`*^9}, {3.51637293901157*^9, 3.516372942537176*^9}}], Cell[BoxData[ RowBox[{"57.1536`", " ", SuperscriptBox[ RowBox[{"Abs", "[", FractionBox[ RowBox[{"0.47`", "\[InvisibleSpace]", "+", RowBox[{"\[ImaginaryI]", " ", "omega"}]}], RowBox[{"41.01`", "\[InvisibleSpace]", "+", RowBox[{"12.8`", " ", "\[ImaginaryI]", " ", "omega"}], "-", SuperscriptBox["omega", "2"]}]], "]"}], "2"]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", "\:0421\:041f\:0415\:041a\:0422\:0420", " ", "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"Plot", "[", RowBox[{"fvihod", ",", RowBox[{"{", RowBox[{"omega", ",", "0", ",", "300"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", "Full"}]}], "]"}]}]], "Input"], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwVlmcgFQ4XxjWUIpmRnYgiSQjJObjXyN6zZCRk3owQInvcZEZmSMnee6RE UcilISMplIyE/JXX++l8+Z3zez6cD88RGzf9KzupqKhO7KCi+v+sMt78PV/y XkFDt4w3ui0QoroujeiR6+B5hLPT5GwQxH8I8echvwDxE+vXi6uC4afwE+ro ogHofDn/92laKCR9SerbTR6BGf3nv2PyIkEuJzA5yOMTtBleiFx5FQeLJ6vF AoomQUsgLY+4JwE6aR8nyGV8gZmDszEvSpMhZTZzdS3uG8jHOgdSfUwDxxcJ 5jVBc5Aqxy4/GpYF8gURrSSPeajXtSePDeVCpKIC473TizBO9/M0IScfPMpl Or2LliCpbFfV+9xCMOc942V09Bc057c8Fh4uAgL5pJBkxgpwljxripctBTZn /piluDUQCr1lotRWCVQfuc737/kDZbueUXiLq2FOnW2hNGgDfBPOS7c/qoVm IToDZ4+/kM1gpQE9jVCx5Bd3vWgH7lCo/feurgNMSuyMfsbtxa64UKaSSy+h 94t16cULNHhXn68lwu4VIJfV3r49+9A8j6dk064XTsSYNTwJ2o9/21iPvTJ4 DVlPjZk45WlxVjrVT0vxDTBvGFyLXqfFT/Ax5xNzP2w6aHM7ehzAtiNO8SwL /eCRo+H97iQ9rv+lZh82HoCvI2pvVOfoMUuyY+pO4wA45OVi4mkG1PT2GSnw GgR+4Jt/eYcBX1StLM72DsLoh6y0nT8ZcEPfpPg371vQY8pY9ihiRFs2ASur hrdAV8qR/ZiGCcM+n43o3HoLXeppGpP2TOhIqRP5T3EIbk2zrbM/Z8K4N+It DLeGQC44JV/3KDN2X84VYmwcghUuVr3IYGZM0khPol0cgrL6xL9t48zYlF1w ipGfAk6GTEVr51kwYV54XlKXAgKL8canMliQW/C7QYwfBcZiDu66usGCGe8t 1E/kUiBNiFyWZcqKSsNlE8LPKGDYSWc5XMuKceJjotlTFKC3iqGhZz2EOmrB rA+3KNC9sa+GeP0Qily+J2bCPgwhKZHWAQOH0DWbL7H35DDIS+ylrznFhif0 fJ7y4TCs9YU1/ohjQ45d12+56gyDHeyXYVxhQ/mtt8caLYdhsJxcI2XOjoP7 ZFMb7IcB+ZnPmLez430TEY47LsNQmphaHnjsMFYvND2z9BgGLmousbzYw6ih Hm1BIA1DlHfOkxfLh/HATT7Zq+7DsPpN4PgPUw48yOd29uu1bZ/Z44cMbRy4 zncmY9Ju2/fypICUICeq+jHQ+Fps++Qrc81iOHFz5aTDq+18pSXSvIFLnGhy i3duVXEYOHmbMh6YcGF9xM8ldolhiIwHjhctXBj9hJFJm2/bt+NZ6vej3PjJ RYi6km77/pfeBMlFbuTmH9x0n6AAGusxmBnzILW3E9fWCwqUvqDEBTTzoOdP ghJ/KQWiisYiuyJ5Udk3lzfUmwJ8MjegOJ4Pi/HaON3GEIy4ywW2dvJhT9N5 Qs67ISA/3mzuX+VDF2la52u1Q7DJESz32/IIMoekTM65bfN/o6QUjvOjiGk8 p8ToWyBLa3rqWvKjg1L+T9Hat0Bwo6+yucOP2SHOtf133kLlRIJ4xG9+vBLl ZwBK23znfZH+jqPoKCNa8Dp/EIhRJUdszAWxrlt66fnVAahiGTgQHiOMYrMh M5rMr6HZg/OLVJMwJjQJRvr398Gz11cap+eEMapVYlyE3AfDkRv2xAvHcbSj LpOdtg/WN4+27953AjWPC3xVo+0FhWlv0u1wEazdG67fLPgSVJU61M7UiuDx iGyP2997QDeblndqWgTrn5CPjVT0gLVZzisloij+tdfUOIo9ENrXI7Bz90nU cg72V7bthpc1XCO3gsWQb50rjKuzCwaZrpaIV4jheb9dK653u+CjW8XtiQkx rHmqn6p8uQt+HFcRRzyFseaKqW5bz+FgllvU1r9TmBfxhL9W+TmcUJ7khlZx fJW5g0Q/1QnuIpL9LecksKZSypnbqwPGzqpucRtI4FhGkYyefgdoEczFAp0k cPNVV+rEqQ44cfFWrEKaBJqs9O2j/9EO0+RetbZVCRxwX6WOdGgH82W7jvaK MxgUp3ZW6XMrdP/zWeTrOYPL1jrCZwpb4SxtDG/wxBm8L1oy0+zcCiwCFQGK ByXxqcvgSOKfFnht9Ff2qbMkirHPuvJytIByQ3Jlp5AUmlOzact7NUHl80eT R0EKFxTO3ltSagK+wSaGUGMpdEvnbxVkbIJ/c5OuhDApzKe937y7rBHqucRE nn+Wwl8f9fS8FhtA9FZXXlemNPJkTh2fj6yHQ6prST3MMjjiU0Fy2VsLy90a 6wInZbDr1adZw/Ea6LuQY3FLRQb/UxiHvLoaCNVW5z/rK4PcotJ+bo41sGx0 vyx/TAZN7Wca8gaqoc8OXgY9lkVRLcMar6oqePQ18eTHp7KoG3peWv5uFYQ6 zMRLj8qi8eQYmLhVgZzzXZN5ejlMa5QZDxfd5q9PfTH3lENSWnl5eVElhIZE bEnBOWyzna+3rKmAy7tHbe6ansMF03DhP6kVcC5cvOuHxzm0cJa69tGvApaj 3sfm5Z9Ddr5gYVfFbf6uCAfjfnl0SvlsVDtQDvI5byR/DMmjFwfTI73NMvjV wub0wEkBA6KiCA92lcC6nniaZogCNrSxmft8KIbNabXu1TQFVLfP/J5TUQy7 6f2OafQoYA59H/3Py8Vw6NLo519CgCUWsu2mnU9A5l+OBXEa8Ky7zs2mlCIA 3uRjK2WIV/Obq8ikR6DkeEShpwnxW/pFYy/9R0CoKjHKerH9saEiRmkSj0BN tStUbRyxieUYXeqvQtB1W/ucQa+I6j1HlMxuFIJNm1k20UURb86F6gbEPoQw Kx721BNKyOkcLlX/MR8U7F71OJ5VwqYTdj0C7fmw6nDDX56ghKvq3BVd+fng QBocm7yohOpJdfaPXPPhQlh4gehdJWSU1hSl7M4H+ic/JZ6uKeFoRsrhq3J5 cG+1VWv+mTI2VEQXaL7IBd3/rm21DSij4W/NoJsVubCP6nBFwpgyJiRQFb69 nwt++66zyKwrI8GQ3X3VPRcucgl9DBEh4LcFFsbd3LnAr3THgT2BgHI/D+9a 8c+B4jir28pWRGTUqbL5VZAF5sG7TD5cI+J4bltc5O0soPEqPOFxg4hco22P jK2zwNZyYSDrLhFXVrg773NnAYdIEN9GJxGtsy9aNd/LhMjurJZyYRW8Rf1s fOleBtjtGlvjXlbB9foKW4XWdGBaC35V/U8Fo17sjTn8IB3a5wSzNWhVsWU5 z5UvPB04B11UfAVU8f70vr3l2ukwmPM3achYFfXWCL96J9MAFbglYppU0UEs R7j6YBpw+lo6r4eqoUXLoCBLZCqMeX0uqI1Xw5EmbSszUirkkBzGPTPUMHnE Ub3ZMhUEnD31l6rUME+yxPzD6VQQtYqVmfushjpv5Bs8PqWAArF59yiq471X zWyTsilgzciV2bapjsVhpo+99iSDAP2D4QCaC3jzheAq40oSfN0vzCDPcgH/ 5Oiwd00mgdNuydv1IhfQx8RezqslCUhrGo4V5hcwuXw6dc4rCUI/+Uvl1V/A 53NzIwe/J0Lh44994Z4aSBo6rkGcTgA+65BcpWANZPtn46s4lAD32I97/YvT wH7aPTaGnQkQHeHN5V2ogZvW/U9bchPA+QrjtSsfNDDU95xPulUCiB9RpSGg JurV6diaj9+FxnuVSjvotFCthYl3Yz4eJHTNDrWwa6GYw4Rk+Hg8PNlLNXdD UAurGMMyBQbiIcNb6+6ighYaflv6kVodD0GGM+PjHlqo4+gia3QzHggMPAGt I1r4zd7vD8PBeHgdEVXnn6eNy9N/X9wLJoOFGtPSpTJtvNJx6o+/ExlmaO6f UGrSRnJcg4SvARl2Rhdn0rzVxq05YYEOQTJIx765nbxTB3t3naeEbPfv7HhW 3TJrHfxv8NJ6MWcc1Pfqei4d1EXhFrbvBX0x8E+UlXGKVRddMxbFQhtigBj3 vmSIUxfPn+ziv1EQA4NaNt/qhHTRPupAdk5ADPx4fd08EHQxtMg5ZVEsBo4M pCCtmy425anGk5KjIZoySif4WhcHGOS9ZT2iQKMko3H9rS7u7WZ0+WkVBXRh lg6973XRV/7NVLV2FJAlRzuvT+vipyLXX16iUZCY+NHv6aYulg2VsnTNREKG 3ocZKxE9HPmn9jDaLhJK+kaepUfqYan4kUx6hwhwLUgluZL1sLqvYqXbNALE A0z4lJL0cJRWXSBBPQIqRUf8Z7O39wWU2nRFIqA2ZlhCpk4PqU3/Dt1cCIdW dUoOZVoP/ZXuXrnlFw79zwdvMijr4y6NK4ca08Pgx8aDNRt1fXzSIa+yJzYM aMSvk2p09BFRxMo6IAwU05kdzC31UdtNb1P5chhUuhoa5HvpI40Wb1P9sTBI Zhs+LvNIHzNtWt10akPB3PH9sNUBAzSfGD5VMX4bvtBNiJeNGGDFM2TPMgmB aJXcS8mfDNCIo7R1TDMETt+yifWfMsDRwN89wkohEPjryzfVBQMU1EiZeica AuwfZrMm9hjilz1HVwt3hoBW4S86JmlDtFiNsFkvC4Zaxb2znkmGaNNb8JiT MRgivMVy5fSMsPZkfPTz4SCIPECVkW1ihKURIwlznUEQnT+QsvuSEU4FU4h8 FUFAHvSMfe1khMTqqYstMUGQLNbkYxNqhE3MW7LNikGQ91VVO7rOCPnsejia ywOhzdj6v/fcxnh8OaOWmB4Aa9JJRr5zxsh8v/36h1R/8CRaW7cvGqNNtFT/ t0h/WDQQc9m7Zow//pQRdvj5w4x7d2jyLhNkiB9QNLX0hw9Fm1UVXCbIomMg mHzEH5p5rjDNaZvgvXSxOCj1g+A90m/Mq0zw6L9p1p39vkA7PKImf9MU6cRq VHokbgDTbPJethBT3OPzlcFQ6Aawbxp0LUWYIuOf26tfOW/AMf5+5cJEU7wT mvBbgPoGoGu3AtMTUzxvNLxA9c4HPPc0SM68N8Ubb4xMJYJ8YFQ6nS/prBk+ XIt8rfbWG0pTLde/L5shVwfPp954L+ApX/+uHWWBQjXLHB0nrgP3d8r9rN0X MS8jkMOR6AFzN1/eyc2+hCy5YxNvaNygxFvbOunwZSxeTwhuZHWGk/jO9DXl MtZN3Nx/utIRmGtclidzrbFvwTdD0u0qKDYdnbt22QaLt2jNJpOuQEe52g4Z TlvcOdSuYLtuC4V6vo2fpm0xn16I5p+TDVBuvxz488AOFViubvjwWYOTtVDM gtsV3CWT3bPSZgV6Yr9shWTtMaX5DF1u7EUIih1Pp6zYY1ceJBlWWgBLysj1 i+1X8bnn0/OFE2awKL5YZBHogM8utajxGJsCT5PPjgwZR+SRNgpPXDGGqsX9 GzU7nVAnqZ5zPtkI7tS16qq0OeHXY8wNmk6GMHg8UL085hqmOU8FnUYDsDx5 KSf+gjN+P2AwrYT6IMNXcqiK2gW/fbENlTbVA//WvafK+l1wkJXCHuamC9Q/ QwqOkF2xb7LUaFZTB1a797W6WLghhcmL9229FmT0bG0EMLnj0wcbj7jUNcF6 bb/KLXf37XwGD9L2acD7DU4bt2F3JAfbNvt0qcNvlfMlaVIeSD5sFNtdrAaN +w/RTCd44O8fhCvdj1VBIn/lDPe6B3Ycu3TH9oEKaIdPkeONSbhQ3fDnVRMR VC+QNaVNSZgZsNozXUcEpJfdP2pGQk+qgodU1UQ4nUIOE7pIwptpw7dli4nA 8lD2ZqstCc0UGTvrMojw8fkdxx/uJFSsT9vxNYAIDtTyBLUYEp6tJM7vUSSC dc/XHT9jSdgdfCdJ8TwRzOPutiWSSXjkTop8gCwRNFm/yY3fJeF0pnvc2mki nBJMOO19j4SlS07C6/xEWCXM8OQXkHA8v3RFmJoICzSJo+qFJCT8Cwi4voMI M73n0xcekXBngsuOtr8EeG+YyHqumIT/xWj8MV0lQLOdAt1gJQkdlC9l3/tG gBrh2R6fahKqjEStzU0RoPRHYgR3LQnnR1nVFCYIkOM5u9OxgYSq0sv9394R IE0mqZ2+iYTM9kk0ChQCJGwqBFY3k/BBm9zZ5AECRLfPnjNvJaFzPNlyvo8A t0OT/my1kXBAQ/aGyksC/A/cpa06 "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, PlotRange->{{0, 300}, {0.0006344646539652115, 0.3507243055293158}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"p", "=", RowBox[{"NIntegrate", "[", RowBox[{"fvihod", ",", RowBox[{"{", RowBox[{"omega", ",", "0", ",", "Infinity"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.5160898890711718`*^9, 3.5160899559640894`*^9}, { 3.516343711824143*^9, 3.516343713181345*^9}}], Cell[BoxData["7.051581730451964`"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"omega0", "=", "0"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"shag", "=", "0.1"}], ";"}], "\[IndentingNewLine]", RowBox[{"While", "[", RowBox[{ RowBox[{"shag", "\[GreaterEqual]", RowBox[{"10", "^", RowBox[{"-", "5"}]}]}], ",", RowBox[{"{", RowBox[{ RowBox[{"While", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"NIntegrate", "[", RowBox[{"fvihod", ",", RowBox[{"{", RowBox[{"omega", ",", "0", ",", "omega0"}], "}"}]}], "]"}], "-", RowBox[{"0.95", "*", "p"}]}], ")"}], "<", "0"}], ",", RowBox[{"omega0", "=", RowBox[{"omega0", "+", "shag"}]}]}], "]"}], ",", RowBox[{"omega0", "=", RowBox[{"omega0", "-", "shag"}]}], ",", RowBox[{"shag", "=", RowBox[{"shag", "/", "10"}]}]}], "}"}]}], "]"}], "\[IndentingNewLine]", RowBox[{"omega0", "=", RowBox[{"omega0", "-", RowBox[{"shag", "*", "5"}]}]}]}], "Input", CellChangeTimes->{{3.5163435178690023`*^9, 3.5163436794072857`*^9}, { 3.5163437220577607`*^9, 3.516343761947031*^9}, {3.5163438643456106`*^9, 3.5163438797740383`*^9}, {3.516372786427702*^9, 3.516372865644641*^9}}], Cell[BoxData["161.93358499999513`"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"p95", "=", RowBox[{"NIntegrate", "[", RowBox[{"fvihod", ",", RowBox[{"{", RowBox[{"omega", ",", "0", ",", "omega0"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.5163436706244707`*^9, 3.5163436719036727`*^9}, { 3.5163728835690727`*^9, 3.5163729252991457`*^9}, 3.516373024967721*^9}], Cell[BoxData["6.699002619670775`"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "\:041f\:0415\:0420\:0418\:041e\:0414", " ", "\:041a\:0412\:0410\:041d\:0422\:041e\:0412\:0410\:041d\:0418\:042f"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"Tkvant", "=", RowBox[{"0.1", "*", RowBox[{"Pi", "/", "omega0", " ", RowBox[{"(*", " ", RowBox[{ RowBox[{"\:0433\:0434\:0435", " ", "0.1"}], " ", "-", " ", RowBox[{ "\:043f\:043e\:043f\:0440\:0430\:0432\:043e\:0447\:043d\:044b\:0439", " ", "\:043a\:043e\:044d\:0444\:0444\:0438\:0446\:0438\:0435\:043d\ \:0442"}]}], "*)"}]}]}]}]}]], "Input"], Cell[BoxData["0.001940050085095002`"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ht", "=", RowBox[{"InverseLaplaceTransform", "[", RowBox[{ RowBox[{"Wt", "/", "s"}], ",", "s", ",", "t"}], "]"}]}]], "Input", CellChangeTimes->{{3.5163731554551506`*^9, 3.516373156203952*^9}, { 3.5163731948256216`*^9, 3.5163732020172343`*^9}, 3.516374013039666*^9, 3.516374301290577*^9, 3.5163744371980157`*^9, {3.516374549580613*^9, 3.51637462187114*^9}}], Cell[BoxData[ RowBox[{"7.56`", " ", RowBox[{"(", RowBox[{"0.011460619361131432`", "\[InvisibleSpace]", "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "6.4`"}], "-", RowBox[{"0.22360679774995673`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]], " ", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "0.005730309680565717`"}], "+", RowBox[{"2.0720569441832595`", " ", "\[ImaginaryI]"}]}], ")"}], "-", RowBox[{ RowBox[{"(", RowBox[{"0.005730309680565717`", "\[InvisibleSpace]", "+", RowBox[{"2.0720569441832595`", " ", "\[ImaginaryI]"}]}], ")"}], " ", SuperscriptBox["\[ExponentialE]", RowBox[{ "0.44721359549991346`", " ", "\[ImaginaryI]", " ", "t"}]]}]}], ")"}]}]}], ")"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{"ht", ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "10"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", "Full"}]}], "]"}]], "Input"], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJxF1Xk4VVsUAPDrXrMipCQZGsgYJZmyNkURJSEhSjKFENKkgZIikilCpEmG iIubeYjMIfM8ZZ4umXnn9Pre++t8v+98e5219l57HUHzq7qXiQQCoZeGQMCf AU4P+VKKroA8kfnzSMpLaPVgKBzLvwHxoiVST23eA2sg+TNz8wOYFeOe86j/ DJH5fvx0Aj7QtKytNkgig/ikxbN1K39olT3NKknzFTpYo9ZWmgKBGpV/ra4y D+wOu9rP87+EUIqMqbttMXxp6I6bPfIKSjoUjbc9+QZLV7Tap62igWgvufD5 RRkcrXq+53dTLHgl+ZobWVXCUNixglmbOGji9QlJ2F8NvharJtTlt3DRbY+9 20gN1C9bB03xf4Qv7/ca9wb8ALdSPqnJlHiocJeMlTldBzwvGirGjyRA/DvT iV7GejAXRaRRq2RgXUpTMzNrAOLyzZrZ25/hEe+hvKS1Bmjo1t8xq5YCysvF UudDfoJ7AkvWTFMqVGhMGZcmN4Jm4CD9zOsv0MVpcmxeugl43Qv0pm3SQOxx SKxDYhPkH7k+NbmcDvpv/U+zPG+GQBFd5cliMribiHWELTaDBZuE74RfBvjp vjZSNmkBhrZe4XH+LHA7EnIjcFMrNOfnuI4NZUHWeWWPKPNWLO+wotEUCjwU Tz1xJ7kVbvteYx+9+RVINx+6GM23wknnk2YjR7KBQ2Dt3AWFNhAwFEkc3pAD LkA5kOneBjOHaZeHfuZA6C/hrU9S2yCEiRLyyyoPnqg4/pDd1g5p1zWdr3Hn g+TarLWJeju8Gjm0jXo7H9AD+ZfPr7bDw/N78p178kFEOV52MagdHGo5rGbU CuACNYs3ndwOZ48QWJ3jC+Abu5TMYEM7IPJ42jRrIXxz0LmRMNkOIiJtxk7X CmFWoYNGmKEDOF6VEaebCoHANUt02d4By6zkj45KRXB2c2njZ/EO6L//Rmfq dRGEvt5PXFLogKrZgPmrdMVgOiilcEW9A8hWHlGTNsVAly5mw3+qA6Jbr6hd rS6GLF+JNCn9DnisfW5sYn8JSEfP3M007ACnfPUXDqElkH+/1bXmXAcYHZBR mFguATPNr4/9sPdH3gn22F/4BrFVIYtM2HrxbWyPx4u/gWbyfhNjLD6X74qk vUgpzM+8iAg91gFra8M/x/xKYZ1sulZ/uAN+OTXdtpsphf5cHrE9+zugtr94 15hBGWxwPC8evbsDMs+mll/5WgZFnS0u+lwdEFMe7TTK/x1yHp3tM6LtgNmw dNEu3+8Q+En1N99UO6hbVfTVLX6HygYFG47WdhgjzetT6svhvWOBROeHdoC6 DWxJqAJ67a7zdPu1Q+DrnWUxiRWgS/753sKpHeQOn1Tw8a6EV2nXPVVk2uEp iwX19mwlqBlkS+ZytkNny40Ex4tVIJc0fJh+pg283N7yGSpWwxf1d1kKn9qg NnmFKDxZA1zEYMoAextwAJPpRslaYKwrV3cebgVFoSAFb4ta8H38phTlt8Kz 2Y/Umz9qYZrtZvjPK60g8/znZfNPP6C//dTwruwWuFsuriltVg+XNzG4HNBs hg8pGUKfguuB+75aYx9PM9SFqZL2VNZDferQmuNoE+y2Mszmlm8Al7HtVUq+ TVBO8pIkcP6EICfLBI2qRthyuI2j5lsjuI1fCBo+/ROSkn3a7CRbYde0Vrjt vToY4Rjx1bzYCkSrzodPDOtA2E1TeW9QK9h9vMVKlaqD10osMX2LraAu9FIj tvcHBH73tTQqaYNP30WrXmhg86PPf/qYSQeYC2+XNheqhVT1qdg9/h3gLLat /QqpFiY/6uiRCjvg2ZZP+QEFNWDtyJ6RK9wJRUl3787J1YDRauDtg9ROWO2+ 5h6wrxqiWSo82dW6QVF5OJNRqRKUD98yuWfaDU/5fBTuEyqh3UHs4OT1bsh9 GTMfWVIB3HVPByvju+ESTVxXwekKeB6mpfF4Uw/ckfTXjHEqB8891azrbT2g zG+YvammDHae9fhlP9cDqnYzU66RZZD/WDK/nbUXHg3ftXa2K4OV0WdOFJVe 2OhnuCK+oQxcU081uL7vhdG2tW9PdEvBCn68HHPugwqLkt7xxRIgxP6mtD3t A0oFlxprXQmE0fK2l8f1QTBJIz4pvgTKvlvyfWzsg2MV16dtz5fAXr2V2MsK /cDRDo4pZcUwZC2c0EkzAFxC4+ku6UVwr0KrqopnAFrC7ctYgouAW9J5IvvA ALTpiHptci2C49RsqYjLAzC0QcnGVLYIPt7RTT/7fQCKj1ccGsspBNvAO7m1 AYOwyVN94mdXARDnYrvyPg5C3yF2zoKSAgg/W0ZILhyEuJQ1UcaEAijn5Tzi NzsIMvXxC29vFIDY+w+lGoa/wP6oUAVxWwGMfa2vLeQfggc8b7XOW+WDF9/i dIrcELy4JGjhq5sPvPf5OGNOD8HP4TZ9BuV80FK30b/rOQQCGh70HVz5kFi7 1qL4awi2BMw8aXueBw4DIn1pScNwqlYnP+t6LnSbtjAWfB+GfqdSudkzuXCm +bFkVf8wcMkKjeRJ5YJcxa8bAzwjsLzR6CtxLAdIKW83bfEeAeVkt4xV6xxw E9WX3Rk7AqR9RAcfjRwYekNrIpkzAp+9P94KEc2B6pBL79WpI1C6FsItMZEN YbcFla+bjcKchu0y961sYP5de8nz5ii4tYgN3DLLhjtX7/n4B49CFznh8+mj 2WBu3tXwvmIUhEttt59iywbx45FXmg+Nwbj27S/m8V+hgHPbS/lN4yCY79zq Q6WAjF9ZrprYONwwrtY+10mBd/Tu/afVx8ElKjXnwXcK+C427rO9PQ42NN1M na8pYNAV9C18aBwUY1NfBZyhQJnh0bF3pAkQtOXdzqxCAcU6KvsXvgmwiPFY rZOkgECJ7vkKvQkAsTcL2iwUGInfNLucPwF838roS0uzwGR3/jbG9gnYHJdV 9DIjC2oir8Lm+Qkw7lu3SX2fBWkB1U/EJSbhwht2upTHWXDXzVfw/MtJEKFJ PFuinQXMgl7CXXJTIOJoL2PZlwllSs/3eahNgUPLmZ2jTZnwyDDq0A7dKeC6 2M8dXZkJxIDMY8ZXpsDOPlj9OTkTltbGrJoip0DxlQE89c2ETJ6lq27xUzBi arvs8iAT3GQZ3LkypiDIeVr7sXsmzNgLPj5TOwVlneuiIpczYbRd/0MNcRqs 53tcVVEmxC+Yf3Zgm4bhhYyHVocywXqzY+ZG3mkgFhF0UiQzof/Ek7ITB6dB gMGm5vOOTGin5A6VWk1DRNtWVp3VDKgMExbJr5yGcAUKObwwA56kyUibtUxD z+Sgf8XXDDheqyK/NjANm9Rf/eZOz4BiBhMNpfVpkLUdjmH4kAHZbgE2mdIz cLlI+che/wy4GRjpdFZ5BuiU35yS9snA/gvxN35rzsAebZm7ZzwzIG2g2EfG YgYaWZlTeq5nQKLe4sfPITNA4TkZxGOeAXaO9Kmn3syAxP7evh6TDBD15aRM JM/A1/KYU3lnM+BdkUS5+PcZeNVjyPlVOwOi9puPfFiegfF9GjkJChnwbFOF aMwFKhhED1yW4s6ArEM70k86UMHtExLs4MyAftOrsHKLCvSGWl2RbBmgmMip dzaUCpea5lfPMWTAkIbJnY3VVEj+StL9Mk+GzU5JDF/bqPD0nnPcbyoZIIwQ aD1MBbVHJnonpsgQMhj3roh2Fr7TewQJD5PhiOd4zQ3FWSi9Nqf/po0MDvFg JKwxC1m7kzVtmskQ/uN5f4PBLPze3qut/JMMU/yyi/ucZ+FlhsfVzTVkiMr2 2Dn4YRau6p6asCkiw8Icm6vuljnYblzS65lEhl07zNfXd80Bh56J4I4EMpw8 muaTKD0HoeJf4ws+kuFtoGEUo/YcjPDerxF7SwbdfTGluZ5zkMthNeHxigwe BtTTds/noI5vhck+nIzNR7X2bdFzoMWnOGIRRgZC5fCUC2UOIp3N1ZyCyJBo vZ9HbHoO7pFU4ud9ydDs7xXXvDYHxncsRPY9JQNtRqPkow2/IdhHnueqDxmM 6G4d6RH+Dc9DD6czPCIDY2yRXajpb1h1PZcWepcM5m16eaTK35DXDBs3u5CB OaLH3qnlN1yWnG0rdSZDqpEDb9fgbxCm7sy570QGUqv3DQrNPBiwMh9nuEqG uGbKAadD8zBMdAo0sCWDVtix3s6j87B+eJgkZ0OG2bMNAVq683A0U4nMZ00G tabxcWH7edjKJHRn+TIZBn8KfOiMnYf+VjGmqYtkeBacaKD1eR5W/EKpCxfI IKuvQEfJmYdtJV3W9Ji9G86YBzfPwyV5p9eSpmTYW/+IV4t1AWIbe1hTjchg WzsWkHVjAeyZdCxc9MjAEXADhL0X4NT6ynD2GTJQTtFPBAUtAP3p7HpGzCw1 /Ccckxeghy4mKP40GRKqdOmEBxZgoxRl1+ZTZNDz6/oSNLMAT/Xi+d1OkmFF y86cSFgEIVp2pjZtMmhXPszr4FmEefVW+0QtMkyUZ94I0lkEjtlptmuaWD8+ UdtLNF0ERwtr+lENMihr1jVevbIICxExnpcx+38fPXDi0SLY1B3tu3icDFJl fBM02YvgI3pm0U2dDE4lXuZXhZcgvl5IlukIGfj61el/H1gCi6nYA59UyVBB ZIq/jZagyVa84STmPchv2ufcEvhYF/SEq5ChhRJ8L+7pEmQpb7FXQ2R41HJ2 t1jYEtgMQD4VyHBgYVtZStwS1N44eO8NZr+DUWx5OUvgsO4pwYBZ5fO7qJbJ Jbja7RnUeRirr9pa9cLKErhv44sLxRwxLjo4yLgMTzZLTJ/GPCeaLDEruAyy /K0S5UpYP7/NyGbVWwada32RZYpkYI8oaz2SuQwPTzDTLcqTITfriUd58TL4 /hL3LcJ8pVlr5+kfy1CjcYvkj7mE64eN6cgyFMgypu7FfCOgZcGddwW4tEqi LOXIIJQc8Wp97wqwmb+VlcNcX3UePTq4AiIZutHMmCVYer1fnFyBEwpxDSmH yNDzcGRL0r0VUI93s2fC/CwugSLjtwJhNLEiPbJkUCxyMP36cgXeLapHZGEO JlDflqWuwFFW20A7zBq3lmX6+ldAkob+VutBMqQ4bjjDrbEKTSb3M+plyFAZ znkj0WAVjvPIDqRiHizmiVa1WIX2sawvgZh5tomM2t1dBUORTcJ6mL0K1DwL 0lehKwpxdR7A5s2o1keDolWwFT45X4A5k0uvZrR2FXQPn7v0DvO4jfn2LWOr UHV0pdoRswHH3S+2O9fgBFUmiQWzo9KjlvV9a2Cz58VL6n4yPLH0Ww86vAaG rTrLbZjzKBEn8gzXIL1sjjsJ895LmX2c/mvAMeR17yzmpS/TnDnLaxDj5aWc I43N284FeV2mdaBmiIR9wizJSLjwa8s6+P1sPB6O+ZIJawL7/nVwe5vwwR1z Ja3YUUvrdUgnUA7KYx6U3G+77LoOpT+aeEQxE87JBwR4rsPHa36a2zEfTDzW Tolah5kjnF5rUlj9+hYubD/XwadbQ74Mc+bdK+Fxvevgup6m9RVz/UfnfPkp LB8b5tdJmGnvIvpUOgKy3hxME4JZrv6kbrsAAZXbGfhbYY65dWVE6iwB/XQZ 796MWSqfacH7PAFZROaxbMCcT/uBrusSAaWn3TYjYe726xfwcySgeKjwmNmH 3bfX588O+xDQrktLj+swJ/UvW6AAAup01fpVjvmwSLhzaAgBbRCjtSnCbJLa 6Kf2hoDYtvA2pGEOLz5VHJNNQJKP6g69xCzKNPFjoZCAHLQEXwdiztL27Tr1 nYC2Xt0o7Iu5ubFsafUnAcVI+z67i3nLsIq00SQBSVfrplthDtwoE8Wxkwb9 OihEVcK8U7fuk81eGtSbcZFRDnNKiGNWviQNeqFzT+YA5lr+xAYHRRrUWF5Q I4KZdb8QS4UeDWKjVWTgxhzlVsy905gGdSValnFilvhqLnTjIg1KKDOKZMOs dTRaRdiBBi3N7rvPgPmJAbe7pzcNOvd9WmhBEuvnCPLDVj8axJWXeXsW88cu vRfSQTQoxTGsbwpzmfXzpK7XNGgy17JvGDP9LeYBRQoN2is2EduO+UH0yunZ MRrk/tylrxizjlSmkxiVBuXIyssWYuYruPbcfJEG0aSqBudhpvSO1NbSERGh SO06BfPMnpZTiXxEdDFJouQz5nxy0NX+3UR0X/WASxLmZ8d0/LeLEZHy4xaJ BMyiNqXVPoeISOn4sfz3mM0T0rUtdYhI9EmASTRmqcNO9pEGRGRj62EciXmt StyvwYSIeN+/vhSBOXzyTaWqDRFdX2x9Foq5bn/gCf4HRESUz9EKwKyS5aDR /IWIJq5JqHliZtUUtWGjEBHbsELJfcztrQOP1fOJyJZ5n9Y9zO4rJmVplUQU +VHb4Q7mz3Di2PMBImL34SS4Y/aopbcqGyWiiAc0aW6YtS4WPlqfJqL8liF7 V8xDD+S/2a8RkX22LNUZs8A3YTXNrST07fskcsA8YdBncX8HCUm9KRG0x5z9 K8orcxcJ2cqVMthhNmTaUiwkRULHD0/22mAO0KI9QtIgoUE7/frLmAn13Sj7 JgntLmj2McXsNZwbc/M+CQnG8cacx8xEiCTKPyahz8nj2SaYOSTOFacHkxAj RweNMeY93j+OJX8moctlz3PPYo6PTP5gn0FCft7t8waYJdP8mMRzSchilfkg brkejYoPFSQU+6uEoodZU7HwZOwgCemrJ0yfxlx9+nXyhXESighkPo77jLXH Jv5ZEtKL+xqjg9kkWKEugoYWWfE2mJ7C7DCZqh/MS4t28e1c08L84s0bk8e6 tGieVaP4OGZuyv0c9XO06N2zjOO4X9Wa8dFdoEUCXYE1xzC/W93efd+eFp1I M+tXx5x5NujSrce06InVRhE1zEoOzsXy/rTIayoi7yjmAi+dPQvBtKj5k5Ih 7vKUDb9c3tCifbef+h/Bz5/loa1DLi3q7zvOpYr5ws5LFeIltGh1+PRXFcwD cirioxW0qD2B/xLuicur41YttChEZDADYabJd3W6OEuLGLKkXQGzkIulu54o Hdp5sfPQYcwSCl1noqXokMmZ3CUlzDIEw30jsnTIYSonF7eqr8bg3SN0KFyW QRO3aZy43icTOvTQW8ZeEfNl27f7fpvTIdsvnnK47aT4WFRs6JAYfQUt7pvZ bIWNrnSon1HqtQLmkIaZfUR/OhRo19Ijjzky/AqLdjAdUuW/moY77kL/YGgE HfrR1uqNO3XsZ6TEBzqU1LBbGnc1bRbLuQI61GG+zU8OM93Bu7+SqXQoz9PJ 6hDmDcvzhUuL2PfmyGq4OQoco9QI9GhHSMRu3ALa5vptG+gRKU+2XxY/j8tq RfRC9IjvoKUt7iNiuVGnxemR0yEDbdya07I3X+2nRwU3aqRwG97ZK70f6FF/ 2pnFg5hdglmizxvSo4cnzAJw3zL2uvnBlB69Uda/jvuB4Ko+1YIeZT4sMMUd kDixwceJHo1ZzkriTij5cTPtCT0KY/BpkMHcPxdqwJJDj8jHsvbjjlU5cK2j iB4J7HrJh9vMr9o/uZwe/djVzYy7dTfd9zPN9KjockT/Acw/zjgrvqLSo+JO /1e4n0VvNHRYokcaT52e4tYa/eCCaBjQCPHTTdxlD7oS+1kZUIkRnxHu3BQt AQlRBpRR1sSL+/bqL6V1KQak/PE4K255Dc9zPw4xICMVWgLutO6sQFc1BqT+ Qb9/P35f2YTpci8woGcX5T/jDrWnGTkZyoB2UdvO4NbLekUvGMWAQjv1NHBz 0MntosYxoJ19JMD97JWDSWgqFn94XRS3V0VbdVcVA9r4W5cOt+rW66MpDQyo z89sRRrzujkHo1cbAxoWC6fivrl0XGXvMJZvbXAvbkcR8hdHWkZUwvmsELek 6+laVRZGVLZ5hoJ7NH9sbDMHI7qidOcLbkvDXUJZ/Iyos5vuLe7z3gFhNIqM iNOzwge35oDtnUAnRuRn6XsOt/qdXQlH3RlRu5OQHm5VrvbW3x6MaH9F20nc Cmracsa+jEj6aMhR3KJv91F3vWdE/sOK0rh3/nyt8ziNESkc4t+Lm4eOI2ms gBFRms/z42a5PGud3s6IGJY6WXHTBl/+tm2EEd0v5GHAvVLcuMtjnhFds49b l8I8tjurU52DCfmYBU/jruy/o9d8jAl9KI5uxl28eTJFSZ8JsVQv1+HOPnqB LcacCY1fD6jCnRCnWm5zhwlVDxkV4X5Fu+fh/jAmxGk5kYXbz4IBLacyoUUv xRTcHsXDS0VVTMiAUeYjbofdlem+Q0zo/ffyGNwn+wPE+PiYUaTZmyDcx4md tdyHmNEWZyl/3KoCYm6cOszoHtPOJ7iVlN23s9owo63jjg9xHzQpyWd8wIxc rbfcx703zIxl7QszWkzVuoGbt2FRnzrAjA4cC3PGTdKSHO3QYEHsJa1WuEeK zO+VJbIghswBM9y+trw1Dx9tQO/j6fRxS7I37jhithEVvTE+9Sfeh4E8fxFW lJHgp4WbGeYutlNZkfUYw4k/++mrv5e9hg29Yt2rifvYuNH8Tv1NiKk2TQP3 pmaz27n5mxDXqZt/3FxoQWMkzo629xr/8be3m6VyAtjRtSatP05/XGwqMMeO jJOO/jFr6avo+YscyG4N/tiSzrW7upoDcQUf/uN+xYVyexlO1G2p9MeXnG+l bwjnRNHFin/83a6+IGGVE9U4/Ot9VmLVWhc3I3n0r9ezP41dqNqMHPf86xoO cRZXeS70fuu/Lrlwtl4nigv92vSvvyY9iBAnbUHHNv7r7YaNrrcdtyBZln99 i0ZUp7JjCxJh/teFDrpaOzS2olimf83cfvO4Q9pWZPLX7ionJx5xcyOFvx54 JxgUfZsb7f3rMa+eRdI9biT81zPmMWbWnv+bwC8oJv2EGwn9NU+oQGFhKDfa /dc63vxTA6ncSPCvs612aIkPcyO+vw7au402S38b4v7rgHDL0osCPGjzXxsF na/e8IsHcfz1wtG3i9G52/9z34Cqe6E373+OPSy9verijv+c+NCV3UqC7z97 6186QZr73wfbYz56V/H/5zCRiMSqEwL/+W6+X9iJrP+tvV3/bfluwf+s4lBq zO7zv6N5zXdXzv3v3UXvQ04a7/zPARKFonW5/9uSSvvju+iu/2yeMiPg4/+/ A+o2+BWX/u9cqtDi+vr//gez4eSd "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0.1}, PlotRange->{{0, 10}, {0.025092761034559017`, 0.4586065913117885}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"hf", "[", "t_", "]"}], "=", RowBox[{"Abs", "[", "ht", "]"}]}]], "Input", CellChangeTimes->{{3.516374890534812*^9, 3.5163749181312604`*^9}}], Cell[BoxData[ RowBox[{"7.56`", " ", RowBox[{"Abs", "[", RowBox[{"0.011460619361131432`", "\[InvisibleSpace]", "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "6.4`"}], "-", RowBox[{"0.22360679774995673`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]], " ", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "0.005730309680565717`"}], "+", RowBox[{"2.0720569441832595`", " ", "\[ImaginaryI]"}]}], ")"}], "-", RowBox[{ RowBox[{"(", RowBox[{"0.005730309680565717`", "\[InvisibleSpace]", "+", RowBox[{"2.0720569441832595`", " ", "\[ImaginaryI]"}]}], ")"}], " ", SuperscriptBox["\[ExponentialE]", RowBox[{ "0.44721359549991346`", " ", "\[ImaginaryI]", " ", "t"}]]}]}], ")"}]}]}], "]"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "\:0412\:0440\:0435\:043c\:044f", " ", "\:0440\:0435\:0433\:0443\:043b\:0438\:0440\:043e\:0432\:0430\:043d\:0438\ \:044f", " ", "\:0438", " ", "\:043f\:0435\:0440\:0435\:0440\:0435\:0433\:0443\:043b\:0438\:0440\:043e\ \:0432\:0430\:043d\:0438\:044f"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"k", "=", "10000"}], ";"}], RowBox[{"(*", RowBox[{ RowBox[{"\:043a\:043e\:043b", "-", RowBox[{ "\:0432\:043e", " ", "\:043f\:0440\:043e\:043c\:0435\:0436\:0443\:0442\:043a\:043e\:0432"}]}\ ], ",", " ", RowBox[{ "\:0433\:0434\:0435", " ", "\:043d\:0430\:0445\:043e\:0434\:0438\:043c", " ", "\:0440\:0435\:0448\:0435\:043d\:0438\:0435"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"dt", "=", RowBox[{"10", "/", "k"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"hust", "=", RowBox[{"hf", "[", "10000", "]"}]}], RowBox[{"(*", RowBox[{ "\:0443\:0441\:0442\:0430\:043d\:043e\:0432\:0438\:0432\:0448\:0435\:0435\ \:0441\:044f", " ", "\:0437\:043d\:0430\:0447\:0435\:043d\:0438\:0435"}], "*)"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"hmax", "=", RowBox[{"hf", "[", "0", "]"}]}], RowBox[{"(*", RowBox[{ "\:0437\:0430\:0434\:0430\:0434\:0438\:043c", " ", "\:043d\:0430\:0447\:0430\:043b\:044c\:043d\:043e\:0435", " ", RowBox[{ "\:043c\:0430\:043a\:0441", ".", " ", "\:0437\:043d\:0430\:0447\:0435\:043d\:0438\:0435"}]}], "*)"}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{ "\:0427\:0438\:0441\:043b\:0435\:043d\:043d\:043e\:0435", " ", "\:0440\:0435\:0448\:0435\:043d\:0438\:0435", " ", "\:0443\:0440\:0430\:0432\:043d\:0435\:043d\:0438\:0439", " ", "h", RowBox[{"(", "t", ")"}]}], "=", RowBox[{ RowBox[{"0.95", " ", "\:0438", " ", "h", RowBox[{"(", "t", ")"}]}], "=", "1.05"}]}], ",", " ", RowBox[{ RowBox[{ RowBox[{ "\:043e\:0431\:044a\:0435\:0434\:0438\:043d\:0435\:043d\:043d\:044b\ \:0445", " ", "\:0432", " ", "\:043e\:0434\:043d\:043e"}], " ", "|", RowBox[{ RowBox[{"h", RowBox[{"(", "t", ")"}]}], "-", "1"}], "|", RowBox[{"-", "0.05"}]}], "=", RowBox[{ "0", " ", "\:0438", " ", "\:043e\:043f\:0440\:0435\:0434\:0435\:043b\:0435\:043d\:0438\:0435", " ", "\:043c\:0430\:043a\:0441\:0438\:043c\:0443\:043c\:0430"}]}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", RowBox[{"k", "+", "1"}]}], ",", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Sign", "[", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"hf", "[", RowBox[{"dt", "*", RowBox[{"(", RowBox[{"i", "-", "1"}], ")"}]}], "]"}], "-", "hust"}], "]"}], "-", "0.05"}], "]"}], "\[NotEqual]", RowBox[{"Sign", "[", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"hf", "[", RowBox[{"dt", "*", RowBox[{"(", "i", ")"}]}], "]"}], "-", "hust"}], "]"}], "-", "0.05"}], "]"}]}], ",", RowBox[{"Treg", "=", RowBox[{"dt", "*", RowBox[{"(", RowBox[{"i", "-", "0.5"}], ")"}]}]}]}], "]"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{"hmax", "<", RowBox[{"hf", "[", RowBox[{"dt", "*", "i"}], "]"}]}], ",", RowBox[{"hmax", "=", RowBox[{"hf", "[", RowBox[{"dt", "*", "i"}], "]"}]}]}], "]"}]}], "}"}], ";", RowBox[{"i", "++"}]}]}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"sigmaper", "=", RowBox[{"Round", "[", RowBox[{ RowBox[{"100", "*", RowBox[{ RowBox[{"(", RowBox[{"hmax", "-", "hust"}], ")"}], "/", "hust"}]}], ",", "0.01"}], "]"}]}], ";"}], " ", RowBox[{ "(*", "\:043f\:0435\:0440\:0435\:0440\:0435\:0433\:0443\:043b\:0438\:0440\ \:043e\:0432\:0430\:043d\:0438\:0435", "*)"}], "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{"\"\\"", ",", "Treg", ",", "\"\< \:0441\>\""}], "]"}], " ", "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{"\"\\"", ",", "sigmaper", ",", "\"\< %\>\""}], "]"}], " "}]}]], "Input", CellChangeTimes->{{3.5163749539645233`*^9, 3.516374954245324*^9}, { 3.5163749848681774`*^9, 3.51637529920793*^9}, {3.5163753388475995`*^9, 3.516375545998763*^9}, {3.5164180254405203`*^9, 3.516418028732126*^9}}], Cell[CellGroupData[{ Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"Treg = \"\>", "\[InvisibleSpace]", "0.7165`", "\[InvisibleSpace]", "\<\" \:0441\"\>"}], SequenceForm["Treg = ", 0.7165, " \:0441"], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"Sigma = \"\>", "\[InvisibleSpace]", "429.31`", "\[InvisibleSpace]", "\<\" %\"\>"}], SequenceForm["Sigma = ", 429.31, " %"], Editable->False]], "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"\:041a\:043e\:043b", "-", RowBox[{ "\:0432\:043e", " ", "\:043f\:0435\:0440\:0438\:043e\:0434\:043e\:0432", " ", "\:043a\:0432\:0430\:043d\:0442\:043e\:0432\:0430\:043d\:0438\:044f\ ", " ", "\:0437\:0430", " ", "\:0432\:0440\:0435\:043c\:044f", " ", "\:0440\:0435\:0433\:0443\:043b\:0438\:0440\:043e\:0432\:0430\:043d\:0438\ \:044f"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"Nkvant", "=", RowBox[{"Ceiling", "[", RowBox[{"Treg", "/", "Tkvant"}], "]"}]}]}]], "Input", CellChangeTimes->{{3.5163756536545525`*^9, 3.5163756544657536`*^9}, { 3.5163757406559052`*^9, 3.5163758443336873`*^9}, {3.516375944875864*^9, 3.516375945031864*^9}, {3.5163759898039427`*^9, 3.5163759960907536`*^9}}], Cell[BoxData["370"], "Output"] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{"Stextrop", "=", "3"}], RowBox[{"(*", RowBox[{ "\:0421\:0442\:0435\:043f\:0435\:043d\:044c", " ", "\:043f\:043e\:043b\:0438\:043d\:043e\:043c\:0430", " ", "\:044d\:043a\:0441\:0442\:0440\:0430\:043f\:043e\:043b\:044f\:0446\:0438\ \:0438"}], "*)"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Masotklon", "=", RowBox[{"{", "}"}]}], RowBox[{"(*", RowBox[{ "\:041c\:0430\:0441\:0441\:0438\:0432", " ", "\:043e\:0442\:043a\:043b\:043e\:043d\:0435\:043d\:0438\:0439"}], "*)"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"hizm", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", "Stextrop", "}"}]}], "]"}]}], RowBox[{"(*", RowBox[{ "\:041c\:0430\:0441\:0441\:0438\:0432", " ", "\:043f\:043e\:0441\:043b\:0435\:0434\:043d\:0438\:0445", " ", "\:0438\:0437\:043c\:0435\:0440\:0435\:043d\:043d\:044b\:0445", " ", RowBox[{"(", "\:0440\:0435\:0430\:043b\:044c\:043d\:044b\:0445", ")"}], " ", "\:0437\:043d\:0430\:0447\:0435\:043d\:0438\:0439", " ", "\:0438\:0441\:043f\:043e\:043b\:044c\:0437\:0443\:0435\:043c\:044b\:0445\ ", " ", "\:043f\:0440\:0438", " ", "\:044d\:043a\:0441\:0442\:0440\:0430\:043f\:043e\:043b\:044f\:0446\:0438\ \:0438"}], "*)"}], ";"}]}], "Input", CellChangeTimes->{{3.5163760323296175`*^9, 3.5163761120613575`*^9}, { 3.5163761944919024`*^9, 3.51637629556448*^9}, {3.516376478911602*^9, 3.516376482858409*^9}, {3.5163771007434945`*^9, 3.5163771481363773`*^9}, 3.516377436658884*^9}], Cell[BoxData[ RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", RowBox[{"Stextrop", "-", "1"}]}], ",", RowBox[{"i", "<=", "Nkvant"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"Ttek", "=", RowBox[{"Tkvant", "*", "i"}]}], ";", "\[IndentingNewLine]", RowBox[{"hTtek", "=", RowBox[{"hf", "[", "Ttek", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"hizm", "[", RowBox[{"[", "1", "]"}], "]"}], "=", RowBox[{"hf", "[", RowBox[{"Ttek", "-", RowBox[{"Tkvant", "*", "2"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"hizm", "[", RowBox[{"[", "2", "]"}], "]"}], "=", RowBox[{"hf", "[", RowBox[{"Ttek", "-", RowBox[{"Tkvant", "*", "1"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"hizm", "[", RowBox[{"[", "3", "]"}], "]"}], "=", RowBox[{"hf", "[", RowBox[{"Ttek", "-", RowBox[{"Tkvant", "*", "0"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "\:042d\:043a\:0441\:0442\:0440\:0430\:043f\:043e\:043b\:044f\:0446\ \:0438\:044f", " ", "\:041c\:0435\:0442\:043e\:0434", " ", "\:0410\:0434\:0430\:043c\:0441\:0430"}], "*)"}], "\[IndentingNewLine]", RowBox[{"hextrop", "=", RowBox[{ RowBox[{"hizm", "[", RowBox[{"[", "3", "]"}], "]"}], "+", RowBox[{"Tkvant", "*", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"23", "*", RowBox[{"hizm", "[", RowBox[{"[", "3", "]"}], "]"}]}], "-", RowBox[{"16", "*", RowBox[{"hizm", "[", RowBox[{"[", "2", "]"}], "]"}]}], "+", RowBox[{"5", "*", RowBox[{"hizm", "[", RowBox[{"[", "1", "]"}], "]"}]}]}], ")"}], "/", "12"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Masotklon", "=", RowBox[{"Append", "[", RowBox[{"Masotklon", ",", RowBox[{"Abs", "[", RowBox[{ RowBox[{"hf", "[", RowBox[{"Ttek", "+", "Tkvant"}], "]"}], "-", "hextrop"}], "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", "}"}]}], "]"}], ";"}]], "Input", CellChangeTimes->{{3.5163763045188956`*^9, 3.5163763836110344`*^9}, { 3.5163764362767267`*^9, 3.5163764572119637`*^9}, {3.516376577332175*^9, 3.5163765815909824`*^9}, {3.516376919034375*^9, 3.5163769345252028`*^9}, { 3.516377230988123*^9, 3.5163772464633503`*^9}, 3.5163773480819287`*^9, { 3.5163774539125147`*^9, 3.5163774570481205`*^9}, {3.516377487109373*^9, 3.516377487624174*^9}, 3.5163775230050364`*^9, {3.5163775575902967`*^9, 3.516377573861125*^9}, 3.516377612876794*^9, {3.516377799344921*^9, 3.5163777997973223`*^9}, {3.516377893756287*^9, 3.516377909403115*^9}, { 3.5163779812568407`*^9, 3.516378015920102*^9}, {3.516378050988963*^9, 3.516378064888588*^9}, {3.5164179455059795`*^9, 3.516417961870408*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "\:041c\:0430\:0442\:0435\:043c\:0430\:0442\:0438\:0447\:0435\:0441\:043a\ \:043e\:0435", " ", "\:043e\:0436\:0438\:0434\:0430\:043d\:0438\:0435", " ", "\:043e\:0442\:043a\:043b\:043e\:043d\:0435\:043d\:0438\:044f", " ", "\:0438\:0441\:0442\:0438\:043d\:043d\:043e\:0433\:043e", " ", "\:0437\:043d\:0430\:0447\:0435\:043d\:0438\:044f", " ", "\:043e\:0442", " ", "\:044d\:043a\:0441\:0442\:0440\:0430\:043f\:043e\:043b\:0438\:0440\ \:043e\:0432\:0430\:043d\:043d\:043e\:0433\:043e", " ", "\:0437\:043d\:0430\:0447\:0435\:043d\:0438\:044f"}], "*)"}], "\[IndentingNewLine]", RowBox[{"Matozh", "=", RowBox[{"Mean", "[", "Masotklon", "]"}], " "}]}]], "Input"], Cell[BoxData["0.0023403763204312175`"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "\:0421\:0440\:0435\:0434\:043d\:0435\:0435", " ", "\:043a\:0432\:0430\:0434\:0440\:0430\:0442\:0438\:0447\:0435\:0441\:043a\ \:043e\:0435", " ", "\:043e\:0442\:043a\:043b\:043e\:043d\:0435\:043d\:0438\:0435", " ", "\:043e\:0442\:043a\:043b\:043e\:043d\:0435\:043d\:0438\:044f"}], "*)"}], "\[IndentingNewLine]", RowBox[{"sigma", "=", RowBox[{"StandardDeviation", "[", "Masotklon", "]"}], " "}]}]], "Input", CellChangeTimes->{{3.51637672255923*^9, 3.516376738954859*^9}, 3.516377323402685*^9}], Cell[BoxData["0.0023392907391910946`"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "\:041f\:041e\:0420\:041e\:0413", " ", "\:0421\:0420\:0410\:0411\:0410\:0422\:042b\:0412\:0410\:041d\:0418\:042f\ "}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"b", "=", "3"}], " ", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{"1", ".."}], "5"}], " ", "-", " ", RowBox[{ "\:0434\:043e\:0432\:0435\:0440\:0438\:0442\:0435\:043b\:044c\:043d\ \:0430\:044f", " ", "\:0433\:0440\:0430\:043d\:0438\:0446\:0430"}]}], "*)"}], ";"}], "\[IndentingNewLine]", RowBox[{"Porog", "=", RowBox[{"Matozh", "+", RowBox[{"b", "*", "sigma", " "}]}]}]}]}]], "Input", CellChangeTimes->{{3.51637675677009*^9, 3.5163768082813807`*^9}, { 3.516376860369872*^9, 3.5163768651454806`*^9}, {3.5163770298881702`*^9, 3.516377030683771*^9}, 3.516377764806461*^9}], Cell[BoxData["0.009358248538004502`"], "Output"] }, Open ]] }, WindowSize->{1202, 705}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, PrintingCopies->1, PrintingPageRange->{Automatic, Automatic}, CellContext->Notebook, CellLabelAutoDelete->True, FrontEndVersion->"7.0 for Microsoft Windows (32-bit) (February 18, 2009)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[567, 22, 1652, 40, 116, "Input"], Cell[2222, 64, 261, 7, 48, "Output"], Cell[2486, 73, 261, 7, 48, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[2784, 85, 573, 13, 52, "Input"], Cell[3360, 100, 333, 8, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3730, 113, 252, 6, 31, "Input"], Cell[3985, 121, 380, 9, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4402, 135, 367, 8, 31, "Input"], Cell[4772, 145, 28, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4837, 150, 250, 5, 31, "Input"], Cell[5090, 157, 380, 9, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[5507, 171, 342, 8, 72, "Input"], Cell[5852, 181, 6617, 114, 241, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[12506, 300, 306, 7, 31, "Input"], Cell[12815, 309, 45, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[12897, 314, 1253, 33, 92, "Input"], Cell[14153, 349, 46, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[14236, 354, 333, 7, 31, "Input"], Cell[14572, 363, 45, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[14654, 368, 628, 16, 52, "Input"], Cell[15285, 386, 48, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[15370, 391, 403, 8, 31, "Input"], Cell[15776, 401, 925, 25, 38, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[16738, 431, 192, 5, 31, "Input"], Cell[16933, 438, 9668, 164, 232, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[26638, 607, 178, 4, 31, "Input"], Cell[26819, 613, 932, 25, 38, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[27788, 643, 4943, 131, 232, "Input"], Cell[CellGroupData[{ Cell[32756, 778, 216, 5, 23, "Print"], Cell[32975, 785, 208, 5, 23, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[33232, 796, 777, 15, 52, "Input"], Cell[34012, 813, 30, 0, 30, "Output"] }, Open ]], Cell[34057, 816, 1563, 36, 72, "Input"], Cell[35623, 854, 3120, 74, 212, "Input"], Cell[CellGroupData[{ Cell[38768, 932, 745, 14, 52, "Input"], Cell[39516, 948, 49, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[39602, 953, 578, 13, 52, "Input"], Cell[40183, 968, 49, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[40269, 973, 886, 23, 72, "Input"], Cell[41158, 998, 48, 0, 30, "Output"] }, Open ]] } ] *) (* End of internal cache information *)